Vacuum engineering at a photon collider? 
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Abstract 



The aim of this paper is twofold: to provide a rather detailed and self-contained 
introduction into the physics of the Disoriented Chiral Condensate (DCC) for the 
£N) | photon (and linear) collider community, and to indicate that such physics can be 

searched and studied at photon colliders. Some side tracks are also occasionally 
followed during the exposition, if they lead to interesting vistas. For gourmets, the 
Baked Alaska recipe is given in the appendix. 



O 1 Introduction 

m 
o 

The twentieth century witnessed tremendous progress in our understanding 
of the fundamental building blocks of matter and their interactions. Not the 
least role in this success was played by continuous advance in accelerator 
technologies. At the beginning of the new century, accelerator-based exper- 
iments are expected to preserve their leading role in the field of high-energy 
physics |T]. 

Over the seven decades since Lowrence's first cyclotron one has observed 
a nearly exponential growth in effective energies of the accelerators by the 
increment factor of about 25 per decade (the Levingston law j2j. By the 
effective energy one means the laboratory energy of particles colliding with 
a proton at rest to reach the same center of mass energy). At that the cost 
per unit effective energy has decreased by about four orders of magnitudes. 
This is indeed a remarkable trend and it was fed by a succession of new 
ideas and technologies [2j: the principle of phase stability, strong focusing, 
high impedance microwave devices, superconducting technologies, storage 
rings and beam cooling. 

However the accelerators were becoming ever bigger and more expensive 
on the whole. We have already entered "the dinosaur era" with monstrous 
machines and the Levingston tendency is slowing its pace. The problem 
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with the circular e + e~ colliders is that the synchrotron radiation severely 
limits maximal attainable energy. It is believed that this technology has 
reached its limits at LEP and no other bigger project of this type will be 
ever realized. Instead the linear e + e~ colliders are considered as a viable 
alternative. Extension of the existing linear accelerator technology towards 
higher accelerating gradients and smaller emittance beams is expected to 
make real a design of the TeV scale linear colliders. Further progress with 
the conventional techniques is problematic unless some radically new idea 
appears. In fact the high gradient efficient acceleration is a tough thing. 
In a free electromagnetic wave the E field is at right angle to the particle 
momentum and no efficient acceleration can be achieved. For efficient accel- 
eration one has to have matter very near or within the beams. Then energy 
considerations combined with the survivability of the accelerating structure 
limits the attainable acceleration gradient (H [2] . 

The proton circular colliders still have some reserve left because, owing to 
the heaviness of the proton, the synchrotron radiation constraint is expected 
only at very high energies. The Large Hadron Collider (LHC) with 7 TeV 
proton beams is under construction now. LHC is a very important high- 
energy physics project and we believe that its results will determine the 
future shape of the field. An analogous collider with the center of mass 
energy about 100 TeV seems also feasible and maybe the Very Large Hadron 
Collider (VLHC) will be the last monstrous dinosaur of this type. 

Other possibilities include the muon colliders first suggested by G. I. Bud- 
ker many years ago [3j. Muons, being about 207 times heavier than elec- 
trons, experience much less radiative energy losses, which are inversely pro- 
portional to the forth power of the particle mass. It seems that the efficient 
multi-TeV muon colliders can be constructed despite the fact that the muon 
is an unstable particle 

But why all the fuss? Are these future very complex and costly accelera- 
tors really necessary? The past research led to the triumph of the Standard 
Model. At that the revolutionary 70 's were followed by decades of the more 
or less routine verification of the Standard Model wisdom - the situation 
eloquently expressed by Bjorken some time ago [I]: "a theorist working 
within the Standard Model feels like an engineer, and one working beyond 
it feels like a crackpot". Since then "crackpots" have developed a string 
theory as the main challenge to the standard paradigm jlj. This "Theory 
of Everything" is full of deep and beautiful mathematical constructs and 
is generally considered as the most promising road towards understanding 
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fundamental physics. The only trouble with it is that it will be extremely 
difficult to check experimentally the predictions of this theory, because the 
most direct predictions refer to the nature of space-time at the Plank scale, 
~ 10 19 GeV, and no experimental method seems to be ever able to access 
such energies in a foreseeable future. So the string theorists are doomed 
to face the fatal question "can there be physics without experiments?" [Sj 
for a long time. Therefore, on the one hand, we have a clear experimental 
and theoretical success up to the electroweak scale, ~ 100 GeV, where the 
Standard Model reigns, and, on the other hand one has a very ambitious 
theory without any clues how to check it experimentally. But what lies 
in between, worth of billions of dollars to spend in future accelerators and 
detectors, to investigate? 

Despite its splendid success, nobody doubts that the Standard Model will 
break down sooner or later. There are several reasons why the Standard 
Model cannot be the final theory and why some new physics beyond the 
Standard Model is expected |6J: 

• SU (3) x SU (2) x U(l) symmetry group defines separate gauge theories 
with three different coupling constants. The conceptual similarity of 
these theories is begging for unification. 

• The family problem - why are there three quark-lepton families? 

• The origin of the quark and lepton masses and mixing angles, as well 
as of the CP violation. 

• Solid experimental evidence of the neutrino oscillations require non- 
vanishing neutrino masses and therefore some extension of the Stan- 
dard Model. However, very minimal extension might be sufficient to 
accommodate neutrino masses. 

• The strong CP problem - why is the allowed CP violating #-term in 
the QCD Lagrangian very small or absent? 

• The hierarchy problem - why is the electroweak scale so different from 
the Plank scale? 

• The cosmological constant problem - why gravity almost does not feel 
the presence of various symmetry-breaking condensates? 

But how far is this expected new physics? The logical structure of the 
Standard Model itself hints that quite interesting and crucial things can 
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happen in the realm of the next generation of the future colliders. One of 
the main guiding principles of the Standard Model, which plays a key role in 
the theory, is gauge symmetry. The historical roots of the gauge invariance 
are reviewed by Jackson and Okun [_7| and the review embraces about two 
centuries. In fact one can go even further through history, another twenty 
centuries or so up to the times of the ancient Greece, and find the roots in the 
most widely known theorem from Euclid's "Elements of Geometry" : The 
sum of the interior angles of a triangle equals 180 degree Euclid deduces 
this theorem from the so-called parallel axiom. All efforts to avoid this 
sophisticated axiom failed and finally led to the discovery of non-Euclidean 
geometry. But we do not follow this track. Instead, we start to generalize 
Euclid's 180 degree theorem step by step |8j. The first step involves the 
concept of exterior angle: the interior angle a and the corresponding exterior 
angle (3 are related by a-\-(3 = n. Then the theorem immediately generalizes 
from triangles to arbitrary polygons: The sum of the exterior angles of a 
polygon equals 2tt. 

Let us now consider a triangle whose edges are not straight lines but 
some plane smooth curves. When the unit tangent vector is transported by 
a length Al along the smooth curve, it turns through an angle Acf). The 
limit of the ratio A0/AZ, when Al — > 0, defines the geodesic curvature of 
the curve. Therefore, for such a curved triangle the 180 degree theorem 
takes the form 

i>— ♦ /— - • 

where the integral is along the triangle edges. This follows from the fact 
that any curved triangle can be approximated by a polygon and then the 
total turning of the tangent along the edges (the integral geodesic curvature 
of the edges) is given by the sum of the corresponding exterior angles. 

One can define the geodesic curvature by using normals instead of tan- 
gents, because the normal rotates exactly as the tangent does when a point 
moves along the curve. The advantage of using normals is that one can 
generalize the concept of curvature to surfaces which have no unique tan- 
gent direction but the direction of the normal is still well defined. The 
corresponding generalization is called the Gaussian curvature jSj and the 
180 degree theorem for a general triangle on a curved surface looks like 

E e X , aIlgl e S+ / g eo, CUrV . + //G_e Ur _ 2 . (1 ) 
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Finally, let D be a domain on the surface whose boundary dD is formed 
by one or more sectionally-smooth curves. We can triangulate D with tri- 
angles which have geodesic inside (not belonging to dD) edges. For each 
triangle we will have (JTJ). If we add these equations up and rearrange the 
angles cleverly we get the Gauss-Bonnet formula JB| 



^ext. an gl es + J geod. curv. + J J Gaussian curv. = 2* X (D), (2) 



dD D 

where x{D) = v — e + /, v being the number of vertices, e - the number 
of edges, and / - the number of triangles in the triangulation; x{D) is the 
topological invariant of D called its Euler characteristic. 

The Gauss-Bonnet formula (ED is indeed a long way from the 180 degree 
theorem, but the potential for generalization is still not exhausted. The 
ideas of Gauss about the curvature and the geometry on the surface was 
further generalized by B. Riemann. It was soon realized that most proper- 
ties of the Riemannian geometry follows from its Levi-Civita parallelism, an 
infinitesimal parallel transport of the tangent vectors. The important con- 
cept of the Levi-Civita connection emerged. All these is the mathematical 
basis of Einstein's general relativity. Further generalization of the concepts 
of Levi-Civita connection and curvature to more general, than Riemannian, 
manifolds lead to the notion of fiber bundles - the mathematical basis of 
the gauge field theories |9J. Even magnetic monopoles are related to the 
generalized Gauss-Bonnet theorem [TP] , 

Therefore, both general relativity and gauge theory can be considered as 
stunning generalizations of the 180 degree theorem of the Euclidean geome- 
try! However, returning to the Standard Model, this is not the whole story. 
Gauge symmetry is important, very important, in the Standard Model. But 
the real shape of the world is determined by its spontaneous violation. Then 
a big question is why and how the SU(2) x U(l) gauge symmetry of the 
Standard Model is broken. So far the phenomenologically adequate answer 
to this question is given by the introduction of the SU (2)-doublet of scalar 
fields, the Higgs doublet, whose couplings and vacuum expectation value 
determine fermion masses and mixings. However there are too many free 
parameters, not fixed by the theory, indicating that in fact we do not un- 
derstand what is going on. That is why the discovery of the Higgs boson 
and investigation of its properties are considered as having the paramount 
importance. 

At this point photon colliders enter the game, because in the 77 collisions 
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the Higgs boson will be produced as a single resonance. The idea of photon 
colliders was proposed many years ago in Novosibirsk |HJ Q2] . You have to 
have a linear e + e _ collider and a powerful laser (several Joules per flash) to 
realize this idea. High-energy photons are produced by Compton backscat- 
tering of the laser light on the high-energy electrons near the interaction 
point. After the scattering, the photons will have almost the same energy 
as the initial electrons and small additional angular spread of the order of 
inverse 7-factor of the initial electron. This additional angular spread does 
not effect much the resulting 77 (or je) luminosity if the conversion point is 
close enough to the interaction point. The 77 luminosity can be made even 
larger than the e + e~ luminosity at the same collider by using the initial 
electron beams with smaller emittances than allowed in the e + e~-mode by 
beam collision effects. 

The detailed development of the photon collider idea (T3J, El [T3j showed 
that their construction is a quite realistic task and requires a small addi- 
tional (~ 10%) investment compared to the linear collider price. The solid 
state laser technologies with required pulse power and duration already ex- 
ist. A free electron laser with variable wave length is also an attractive 
alternative |Tfi] . 

The expected physics at high-energy photon colliders is really exciting 
and very rich. It includes P21 [13 ITSj : 

• Higgs boson physics, both Standard Model and supersymmetric. Es- 
pecially one should mention the unique opportunity to measure its two 
photon width, as well as the possibility to explore CP properties of 
the neutral Higgs boson by controlling the polarizations of the back- 
scattered photons. 

• Search for supersymmetry. In particular, charged sfermions, charginos 
and charged Higgs bosons will be produced at larger rates in 77 col- 
lisions than in e + e~ collisions. The 7e option will enable potential 
discovery of selectrons and neutralinos. The photon collider will also 
be an ideal place to discover and study stoponium bound states. 

• Exploration of the gauge bosons nonlinear interactions. 

• Top quark physics. 

• QCD-probes in a new unexplored regime. 
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• Investigation of the photon structure - its hadronic quantum fluctua- 
tions cannot be completely determined from the first principles because 
the large distance effects contribute significantly. Therefore various 
phenomenological models need experimental input for refinements. 

• Search for the low-scale quantum gravity, space-time noncommutativ- 
ity [19| and extra dimensions. 

The last item is exotic enough but one should not forget that J2H| "Ev- 
ery time we introduce a new tool, it always leads to new and unexpected 
discoveries, because Nature's imagination is richer than ours" . 

In this paper we would like to indicate that the physical program of 
the photon collider can further be enriched if it is considered as a tool to 
perturb the QCD vacuum. An interesting phenomenon of the Disoriented 
Chiral Condensate formation was discussed earlier in the context of hadron- 
hadron and heavy ion collisions. We believe that photon colliders are also 
eligible devices to perform such kind of research. 

The paper is organized as follows. We begin with the discussion of the lin- 
ear sigma model, which is used as a QCD substitute in the majority of DCC 
studies. The idea of the Disoriented Chiral Condensate is explained and in- 
vestigated in the third section. The fourth section considers the possibility 
of the DCC production at photon colliders. The Baked Alaska scenario is 
examined in some details. Quantum state of DCC is explored in the next 
section. It is mentioned that at photon colliders a direct production of this 
state might be possible. In the last section we provide some concluding 
remarks. The references on the subject are very numerous and we list only 
a few of them. We hope that an interested reader can find independently 
other important contributions which missed our attention. 

2 Linear sigma model 

The Lagrangian of quantum chromodynamics (QCD) looks "deceptively 
simple" [21] • Indeed, it encodes the description of a surprisingly wide range 
of natural phenomena, from nuclear physics to cosmology, and nevertheless 
is given by the very compact expression 




1 




7 



where 



D = d + igA, G^y = dpAv - dyA^ - g[A^, Ay], A M = AJ-^ 

and A a , a = 1, . . . , 8 are SU(3) Gell-Mann matrices. The theory (QCD) 
which is defined by this Lagrangian "embodies deep and beautiful princi- 
ples" and is one of "our most perfect physical theories" [221- However, if 
you are interested in applying this "most perfect physical theory" to under- 
stand the low-energy experimental data, you will not be particularly happy 
by discovering at least three reasons [21] for your grievance: 

• The Lagrangian (JBj) describes quark and gluon degrees of freedom, while 
"correct" degrees of freedom for low energy phenomena are their bound 
states - various colorless hadrons. 

• unlike quantum electrodynamics, gluons have self-interactions which 
render QCD in a nonlinear theory with the corresponding increase in 
the computational complexity. 

• at low energies the effective coupling constant is large and usual per- 
turbative methods are not applicable. 

However, things are not so bad as they look. It turns out that many impor- 
tant features of the low-energy dynamics are governed by symmetries of the 
QCD Lagrangian and their breaking patterns. For light quark flavors the 
QCD Lagrangian possesses (approximate) Ur(3) x C/jt (3) chiral symmetry. 
The corresponding transformations are 

q R ^e l z ° R q R , q L -> e* * ° L q L , 



; A a aa ■ Xa aa 



q R ^e^ ° R q R , q L -> e l * " L q L , (4) 

where Ao = J\. Fates of these symmetries are different. The first line 

corresponds to the ^V(l) x C/a(1) transformations with 9y,A = \ {®\ ± 0° R ). 
The singlet vector current, generated by UviX) transformations, remains 
conserved in the low-energy limit and the corresponding conserved charge 
is identified with the baryon number. On the contrary, Ua{X) symme- 
try is broken due to quantum anomaly. As a result, rf meson becomes 
much heavier compared to other pseudoscalars. Non-Abelian symmetries 
SUr(3) x SUl(3), as well as ?7a(1), are further broken spontaneously due 
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to a nonvanishing expectation value of the quark- antiquark condensate: 
< qRQL 0. Eight pseudoscalar mesons (n, K, rf) are Goldstone bosons as- 
sociated with this symmetry breaking pattern SUr(3) x SUl(3) — > SUv{3). 
In fact these Goldstone bosons acquire small masses because quark mass 
terms in the QCD Lagrangian break explicitly the Ur(3) x Ul(3) chiral 
symmetry. 

Having in mind this picture of QCD symmetries and their breaking, one 
can try to model it by some effective low-energy theory for mesons, which 
are excitations on the quark- antiquark condensate ground state |23l El|- 
One has two kinds of excitations, scalar and pseudoscalar mesons, because 

QrQl ~qq + q^Q- 

Therefore, for three light quark flavors, one needs a complex 3x3 matrix 
field &ab ~ qmQLa to parametrize the scalar (S) and pseudoscalar (P) meson 
nonets: 

<D = S + iP EE ^ {a a + l7T a ) + y (a + 17T ) . (5) 

The imaginary unit is introduced to make the pseudoscalar matrix P Her- 
mitian: iP corresponds to (775(7, but ((775(7) + = —(775(7- 

The effective Lagrangian for the field $ should have the form [25J 

C = Sp («9 M $ + cP$) - V($, $+) + C SB , (6) 

where Csb describes symmetry breaking effects and V(<$>,<& + ) stands for 
self-interactions of the meson field. If we want the theory to be renormaliz- 
able (although for effective theories this requirement is not obvious), quartic 
couplings are at most allowed in V(<$>, $ + ). The chiral transformations (J2J) 
read in terms of the $ field 



U V (1) : 


$ _ 


-> e 2 


vqq 1 2 "v - $ ) 


U A (1) : 




-> e 2 




SU V (3) : 


$ _ 


j Aa 

-> e 2 




SU A (3) : 


$ - 


A Aa 

-> e % 2 


/la ■ Aa aa 
^C^g 1 2 



Therefore, ($ + $) and Sp(<l> +( 1>) 2 are invariant under these transforma- 
tions and the most general form of V(<£>, $ + ) is 

V ($, $ + ) = m 2 Sp ($+$) + A Sp ($+$) 2 + A' [Sp ($+$)] 2 . (8) 

The symmetry breaking part of the effective Lagrangian has the form 

C SB = Sp H(& + $ + ) + c [L»et ($) + L»e^($ + )] . (9) 
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Here the first term describes explicit symmetry breaking due to nonzero 
quark masses. The matrix H represents the constant nine-component ex- 
ternal field: H = 4f /i a + 4r/io- In practice isospin symmetry and PCAC 
are good approximations because u and d quark masses are very small. To 
preserve these symmetries, the most general possibility is to have only two 
nonzero constants ho and h$ [26J. ho gives a common shift to pseudoscalar 
(and scalar) masses, while h% breaks the SUv{3) unitary symmetry down 
to isospin 5t/y(2) and generates the mass differences between 7r, K and 
7/, as well as between their parity partners (the phenomenological situa- 
tion in the scalar nonet is not completely clear yet |25j). The determinant 
term is invariant under SUv(3) X SUa(3) transformations from (jZj), because 
Det(AB) = Det{A)Det{B) and Det(e i ^° a ) = 1. However it violates U A {1) 
symmetry down to Za(3), because Det(e lX ° dA ) = 1 only then Xo9° A = ^n, n 
being an integer. This explicit breaking of Ua(X) removes the mass de- 
generacy between rf and tt [271 EHJ and, therefore, is very important for 
describing the pseudoscalar nonet. Another interesting property of the de- 
terminant term is that it gives equal and opposite sign contributions to the 
masses of the corresponding scalars and pseudoscalars [2E]- Therefore, the 
large splitting between scalars and pseudoscalars is expected solely from 
the fact that rf is much heavier than tt |28| . This is exactly the situation 
observed in experiment. Physics behind the determinant term is related 
to the U A (1) quantum anomaly, mentioned above, caused by nonpertur- 
bative effects in the QCD vacuum due to instantons [29|. Note that the 
% [Det($) - Det(§ + )} term is not allowed as it violates P and CP [30j. In- 
deed, under charge conjugation $ — ► <l> r , which does not change the deter- 
minant. While under parity $ — > $ + and Det(&) — Det{$ + ) changes the 
sign. 

The linear sigma model, as defined by (JUD, (JE| and (EI), has six free param- 
eters to be fixed from experiment: m 2 , A, A', c, ho and hg. Five parameters 
can be fixed by using experimental information from the pseudoscalar sec- 
tor alone, for example [23 EI], pion and kaon masses, the average squared 
mass of the rj and rf mesons 0.5 (m 2 + m^,), and two decay constants f w 
and Jk- To fix the A' coupling constant, which violates the OZI rule |32j, 
some experimental information from the scalar sector is required, for exam- 
ple [25j, the sigma meson mass. The other scalar masses, the scalar and 
pseudoscalar mixing angles, and the difference m 2 , — m 2 are then predicted 
quite reasonably [23 EH E2] • 
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To summarize, the linear sigma model is an attractive effective theory 
candidate for description of the low energy QCD dynamics. Phenomenolog- 
ically, it is quite successful and explains various puzzles concerning scalar 
and pseudoscalar mesons [32j: 

• why the pion and kaon are light 

• why the rf is so heavy 

• why the scalar mesons are much heavier than pseudoscalars 

• why the sigma meson is so light compared to other scalars 

• the pseudoscalar and scalar mixing angles 

• the accidental degeneracy of the ao(980) and _/o(980) mesons 

• the strong coupling of ,/o(980) to KK 

• two photon widths of ao(980) and /o(980) mesons 

In the next sections we will be interested in some qualitative features of 
the dynamics described by the linear sigma model. At that we will make 
further simplification by neglecting the effects of the strange quark. In the 
two flavor case, one can assume that the field $ in the Lagrangian is the 
2x2 complex matrix. However, SU (2) has a unique property among SU (N) 
groups, its fundamental representation being equivalent to its complex con- 
jugate. Owing to this property, two linear combinations $ + t^*T2 and 
$ — T2$*T2 both transform irreducibly under the SUr(2) x SUl^) group 
|24j . Each of them has only two independent complex matrix elements. 
Therefore, it is possible to construct two flavor linear sigma models by us- 
ing only four lightest mass eigenstates 7r , tt° and a. Hence we take 



and, therefore, the Lagrangian takes the form (up to the irrelevant constant 




1 i -+ -+ 



Ti being the Pauli matrices. Then 




term) 




^(a 2 + n 2 -v 2 ) 2 + Ha, (10) 
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here 

A S = A+-A, v = — — , H = h . 

Z A s 

This is the classic linear sigma model of Gell-Mann and Levy jHBj- Its free 
parameters A s , v and H (the strength of the symmetry preserving term, the 
location of its minimum and the strength of the symmetry-breaking term) 
can be fixed by using pion and sigma masses and PCAC as follows [34] . In 
the chiral limit (then H = 0) the linear sigma model potential 

V S = — (<T + 7T - V ) 

has a famous "Mexican hat" shape. Therefore, the chiral symmetry is spon- 
taneously broken because the sigma field develops a nonzero vacuum expec- 
tation value < a >= v (the pion field, being pseudoscalar, cannot acquire a 
nonzero vacuum expectation value without violating parity). The symmetry 
breaking term Vsb = —Ha tilts the Mexican hat and now < a >= a® ^ v. 
Shifting the sigma field by its vacuum expectation value, a = cr + a\ and 
isolating quadratic terms ^f-a' 2 and ^ tt 2 in the potential Vs + Vsb, we get 



^a' 2 and ^ - 

meson masses 



ml = \ s (a 2 -v 2 ), ml = X s (3a 2 -v 2 ). (11] 
The vacuum expectation value gq is determined from the condition 

dV(a, tt) 



da 

which gives 



= 0, 

7T = 



H = \ s a Q {a 2 - v 2 ) = a ml. (12) 

Besides (II ip and (JEJ), we need one more relation to determine four quanti- 
ties A s , i>, H and a$. This relation is given by PCAC: 



d»J 6li = Umin. (13) 

Indeed, the axial-vector current is nothing but the Noether current 
associated with the SUa(2) transformations 

$ -> S 0i ^S 6i . (14) 

In terms of the a and tt fields, the infinitesimal form of EqUH reads 

5a = —TTiOii bixi = a6i. (15) 
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— * 

The divergence of is given by the Gell-Mann-Levy equation |27| 

( x ) - -^L 

where 5C is the variation of the Lagrangian under ( 1151) with space-time 
dependent parameters 9i(x), which equals 

5C = adf.TTi d% - Tr.d^a d% - #tt^. 

Therefore, 

and comparing with PCAC EqJUUwe get 

H = f n ml. (16) 

Now from (1TTT) . (JT2J) and ( TTfij) it is easy to get 

t \ m 2 — m 2 2 m a ~ 3m 2 2 tt £2 
°q = U, A s = — — - — , v = — 5-/^, H = f n m n . 17) 

The precise values of these parameters are largely immaterial having in mind 
idealized nature of the model. In any case, they can be estimated from (1T7|) 
if needed. For example, for m a = 600 MeV one gets: A s ~ 20, v ~ 90 MeV 
and H ~ (120 MeV) 3 . 



3 Disoriented chiral condensate 



The linear sigma model potential in the limit H — > has a degenerate 
minimum at a 2 + if 2 = v 2 (in this limit = and v = f w ). The vacuum 
state, we believe our world is based on, points in the cr-direction, < o >= 
/tt, < 7? >= 0, and, therefore, spontaneously violates the chiral symmetry. 
The natural question is whether one can change the vacuum state by some 
perturbation. The following analogy is helpful here: SU (2) x SU (2) is locally 
isomorphic to 0(4); therefore, the order parameter < $ > of the linear 
sigma model can be considered as some analog of spontaneous magnetization 
of the 0(4) Heisenberg ferromagnetic. Then, changing the vacuum state in 
the relativistic field theory, which assumes an infinite system, is analogous 
to rotating all spins in the infinite magnet simultaneously and is clearly 
impossible. Our universe, although not infinite, is quite large and hence 
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at first sight we have no means to alter its vacuum state: only one QCD 
vacuum state is realized in our world, all other chirally equivalent vacuum 
states being unreachable and thus unphysical. However experience with 
real magnets suggests that this simple argument (as well as virtually all 
other no-go theorems) may point not so much to the real impossibility but 
to the need of more elaborate imagination. In the case of ferromagnet 
it is relatively simple to change the magnetization in some large enough 
volume. All what is needed is to apply an external magnetic field. Even 
such a comparatively weak field as Earth's magnetic field can do the job. 
We are tempting here to indicate one interesting application of this effect 
|35j . Above the Curie temperature the rotational invariance is restored in 
the ferromagnet and there is no spontaneous magnetization - all record of 
the previous magnetization is lost. As lava from a volcano cools below the 
Curie temperature the Earth's magnetic field aligns the magnetization of 
the ferromagnetic grains. By studying such solidified lavas (basalt rocks), 
geophysicists have reconstructed a history of the Earth's magnetic field 
with a striking result that the Earth's magnetic field has flip-flopped many 
times, once in every half million years, on the average. But this is not 
the most interesting part of the story. Investigation of the ocean floor 
magnetization revealed a surprising strip structure. Successive strips of 
normally and reversely magnetized rock lied symmetrically on both sides 
of the volcanic mid-Atlantic ridge. The explanation of this enigma comes 
from plate tectonics. On each side of the ridge the tectonic plates are pulled 
away, one of it towards Europa and the other towards America. Lava, 
emerges from the middle, solidifies, sticks to the plates and is also pulled 
away with the magnetic field orientation recorded in it. So the oceanic floor 
seems to be a gigantic tape recorder for reversals of the Earth's magnetic 
field! This discovery was crucial in recognition of Alfred Wegener's theory of 
continental drift - the idea which initially was met with enormous resistance 
from geophysicists. 

Long ago Lee and Wick argued |36j that an analogous domain formation 
phenomenon is also possible in the case of quantum field theory with degen- 
erate vacuum and in principle there should exist a possibility of flipping the 
ordinary vacuum in a limited domain of space to an abnormal one. "The 
experimental method to alter the properties of the vacuum may be called 
vacuum engineering" [37] . It seems that a new generation of the very high 
energy heavy ion and hadron colliders may provide a practical tool for such 
vacuum engineering. The scientific significance of this possibility can hardly 
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be overestimated, because "if indeed we are able to alter the vacuum, then 
we may encounter some new phenomena, totally unexpected" pH] . 

Disoriented chiral condensate formation is one of the new phenomena 
which may happen in very high energy collisions [38, 39j. In such a col- 
lision there is some probability that a high multiplicity final state will be 
produced with high entropy. Collision debris form a hot shell expanding in 
all directions nearly at the velocity of light. This shell effectively shields 
the inner region up to hadronization time and then it breaks up into indi- 
vidual hadrons. The hadronization time can be quite large [40J and during 
all this time the inner region has no idea about the chiral orientation of 
the normal, outside vacuum. Therefore, if the inner vacuum is perturbed 
enough in first instants of the collision to forget its orientation, then almost 
certainly it will relax back in the ground state other than the cr-direction. 
Of course, the explicit symmetry breaking (~ H) term lifts the vacuum de- 
generacy. However, the corresponding tilting of the "Mexican hat" is small 
and will not effect the initial stage evolution significantly jUj. Therefore, 
it is not unlikely that some high energy collisions can lead to the formation 
of relatively large space domains where the chiral condensate is temporarily 
disoriented. At later times such Disoriented Chiral Condensate will relax 
back to the normal vacuum by emitting coherent burst of pion radiation. 

But how can the initial vacuum be excited? A short time after the 
collision of the order of 0.3—0.8 fm/c the energy density in the interior of the 
collision region drops enough to make meaningful the introduction of a and 
7r collective modes jl2] . After this time the classical dynamics of the system 
is reasonably well described by the linear sigma model. However, initially 
the o and tt fields are surrounded by a thermal bath. So we need the sigma 
model at finite temperature. To reveal a simple physical picture behind the 
phenomenon, we will use the following simplified approach (221 [44j . Let us 
decompose fields into the slowly varying classical part (the condensate) and 
high frequency thermal fluctuations 

<j>(x) = <j>d{x) + 84>{x). 

By definition the thermal average < (j) >th= 4>d and < 5(f) >th= 0. Therefore 
the thermal averaged symmetric potential, which determines evolution of 
4> c i at initial times, until the effects of the explicit symmetry breaking term 
become significant, has the form (we have suppressed isospin indices for a 
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moment) 

< V S >th= j(<^/+ < • 

To calculate < (5(j)) 2 > t h, let us decompose 6<fi(x) into the annihilation and 
creation operators 



Wx) = 1 rai (<" )e ~ ihx + a+{%)eikx ) • (18) 



dk 1 

with cjfc = \Jk 2 + m 2 and (our normalization corresponds to (27r)~ 3 particles 
per unit volume) 

[a(£),a + (V)] =5{k-k'). 

At the thermal equilibrium the thermal bath is homogeneous over the (large) 
spatial volume V . Therefore, the thermal fluctuations are the same at every 
point inside V and < (5(f)) 2 > t h can be replaced by its spatial average 

< (H) 2 >th^ i J dx< (<50) 2 > th ■ 
Substituting here ( TTHT) we get 

< (<50) 2 > th ^ i y* ^ < a(fc)o + (fc) + a + (£)a(£) > tfc . (19) 

We assumed that the chemical potential of the field (f) is small, so that the 
probability of finding its two quanta simultaneously in a unit volume is 
negligible, and hence 

< a(k)a(-k) > t h~ 0, < a + (k)a + (-k) > t h~ 0. 

However, 

< aa + + a + a >th= 2 < a + a > t h a> + ] 

and the second term gives a temperature independent constant. Actually 
this contribution in ( 1191) is infinite and should be cured by renormalization 
(that is subtracted). The nontrivial finite part is 

< >th= t*7 / — < a + (k)a(k) > th . 

V J w fc 

— * — > 

However, a + (k)a(k) is the number density operator (in momentum space). 
Hence, its thermal average is given by the Bose-Einstein distribution 

- N 
< a + (k)a(k) > th = eLUk/T _ p 
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where N = V/{2tt) 3 is the total number of 0-particles in the volume V. 
Finally, 

< (<50) 2 > th = [ 1 (20) 

J (2tt) 3 u k {e^/T - 1) 

In a high temperature limit T ^> m, ( 1201 ) is simplified to 



, ,o f dk 1 f x dx /ni\ 

< W) >(ft= y _ = _ y __ = _. (2i) 

Let us now restore the isotopic content of our theory. Each isotopic mode 
gives a contribution (J2DJ) to the effective thermal potential. However a- 
meson is too heavy. Therefore, we assume tOk/T ^> 1 for it and neglect 
its contribution. There remain three pionic modes. Pions, on the contrary, 
are light and we neglect their masses. Then the thermal effective potential 
takes the form 




<V s > th =^\a l + Tt l + — -v z ) . (22) 
The minimum energy configuration corresponds to 




J>2 

< 17 >= \/ V z — . 

Therefore, the cr-condensate completely melts down at T c = 2v w 180 MeV. 
Above this phase transition point the vacuum configuration corresponds to 
< a 0. In fact, the cr-condensate never melts completely down (for 
temperatures for which the linear sigma model still makes sense), because 
of the ~ H term. However, near the critical temperature this residual value 
of the cr-condensate (which minimizes V ~ ^cr 4 — Ha) is quite small 

/4fA 1/3 
< a >« \ -j-J ~ 3 MeV <C f n . 

Temperatures of the order T c can be reached in very high energy col- 
lisions. Then, in some small volume, chiral condensate is melted and all 
information about the "correct" orientation of the chiral order parameter 
is lost. What happens when this volume cools down? Again an analogy 
with magnets is helpful. If a magnet is heated above the Curie tempera- 
ture and then slowly cooled, it loses its spontaneous magnetization. This 
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happens because many small domains are formed with magnetization direc- 
tion changing at random from domain to domain, so that there is no net 
magnetization. Therefore, if we want to have a large DCC domain, slow 
cooling in thermal equilibrium is not the best choice. Indeed, it was argued 
that in such circumstances the size of DCC domains remains small. 
Hopefully, the interior of the collision fireball is cooled very rapidly due to 
fireball expansion. Rajagopal and Wilczek found [34] US] that in such an out 
of equilibrium process larger DCC domains can be formed. This is analo- 
gous to the quenching technique in magnet production from a melted alloy. 
The physical mechanism which operates here is the following |41j. After 
a quench, the temperature suddenly drops to zero and, therefore, the dy- 
namics will be governed by the zero temperature Lagrangian. If the cooling 
process is very rapid, the field configuration does not have time to follow 
the sudden change in the environment. Therefore, immediately after the 
quench fields do not have vacuum expectation values. Hence, the system 
finds itself in a strongly out of equilibrium situation, namely near the top 
of the "Mexican hat" . The vacuum expectation values will begin to develop 
while the system is rolling down towards the valley of the symmetric poten- 
tial, but this will take some time. Meanwhile the Goldstone modes (pions) 
will be tachyonic: 

ml = X s (<a > 2 -v 2 ) < 0, (23) 



if < a > is small. Therefore, the oscillation frequencies ut = \jk 2j r m\ 
will be imaginary for long enough wavelengths and they will grow with time 
exponentially. The zero mode is the one which is amplified most effectively. 
As a result, a large sized correlated region will be formed with nearly a 
uniform field. When the fields approach the bottom of the potential and 
< o > gets close to its zero temperature value, this mechanism ceases to 
operate. Therefore, a natural question is how fast the rolling down takes 
place and whether the zero mode has enough time to be significantly am- 
plified. To answer this question, one should consider the evolution of the o 
and 7? fields, according to the linear sigma model. During this evolution we 
have 

3% = 0, ^4 = #tt. (24) 

Derivation of the second equation (PCAC) was given earlier. At that the 
axial-vector current is 

n d{SC) 

5fi d(d»6i) M M 
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The conserved vector current J M is the Noether current associated with the 
SUyi^) transformations from (jjj) and a similar procedure will give [27] 

J jj, = 7T X d^TT. 

To make the problem analytically tractable, we idealize the initial condi- 
tions and assume that the whole collision energy is initially localized in the 
infmitesimaHy thin pancake to an infinite transverse extent |46j. Then the 
fields can depend only on the longitudinal coordinate x. Besides, such ini- 
tial conditions are invariant under the longitudinal boosts. Therefore, in 
fact the fields can only depend on the proper time r = \/t 2 — x 2 . Then 
d fi = —4- with the (Minkowskian) 2-vector x^ 1 = (t,x). Therefore, 



7* % U, -» -1* 7* Xn 

JfjL = — TT X 7T, J 5/i = 

r r 



TTG — GIT 



where the dot denotes differentiation with respect to r. Using d^x^ = 2, we 
get 

2 _> ^ d (tt x 7?\ Id/, ^ 

J» = - 7T X 7T + T— = — [TIT X 7T 

r dr \ r I r dr \ 

and ^ 

t^Jty = -— T [IT (J - OTT 

Therefore, the conservation of the vector current and PCAC ( 12^1) imply 

_> T 

7T X 7T = — , 7T 0" — (77? = 1 / t'tt(t') dr' , (25) 

T T T J 

TO 

— > 

with a and 6 as integration constants. Initially, far from the valley of the 
symmetric potential, the symmetry breaking term Htt plays an insignificant 

_____ — * 

role and can be neglected. Then (1251) shows that a ■ b = and the triad 
a, 6, c= axb forms a convenient axis for decomposition of isovectors. The 
first equation of ( 1231 ) indicates that 7T a = and, hence, ( 1231 ) is equivalent to 
the system 

a . . 6 . . . 
r r 

Because of the last equation, the motion in the (7r&, 7r c , cr)-space is planar 

— = k = const. 
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Then, the first two equations give 

k~- 

To simplify the discussion, we assume b ^> a. Then ir c « and the motion 
plane coincides with the (7r&, a) plane. Let us introduce the polar coordinates 
in this plane 

7r& = /sin#, cr = f cos 0. (27) 

Then, (ESD gives 

fe = --• (28) 

T 

The equation for the radial coordinate / can be derived from the equation 
of motion 

□tT = -A s (a 2 + 7? 2 - v 2 ) 7T 

with □ = as the d'Alembertian. In our case, this equation is equivalent 
to 

1 d 



(^)=-M/w K 



r c?r 

Substituting here (127)) and using (1281). we get the radial equation 

/ + 7 = 7& - M/ 2 - « 2 ) /• (29) 

At later times the difference g = (f—v)/v is expected to be small. Therefore, 

(f 2 -v 2 )f* 2v 3 g 

and (1201) reduces to 

+ - = "xt - 2X *v g- 

r -u 4 r z 

Introducing a new dimensionless variable s = y/2X s vr (note that in the 
H = limit Tfifj = y/2X s v), we get the inhomogeneous Bessel equation 

ds z as \ v J 

Therefore, the solution can be expressed through the Bessel functions Jq{s) 
and Yq(s). Hence, for large proper times it will exhibit a damped oscillatory 
behavior. For example, for large s ^> 1, 



MS} „ yjl COS (, - l) 
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Large enough compared to what number? The inhomogeneous term in ( 1301) . 
which is a reminiscence of the influence of the angular motion on the radial 
motion, is characterized by a dimensionless number b/v 2 , which we assume 
to be much greater than one. Therefore, the asymptotic value of s can be 
estimated to be s ~ b/v 2 , which translates into the proper time 

T « ~ viy- (31) 

This gives us an estimate of the rolling-down time. 

Let us now consider the process of formation and growth of correlated 
domains in a scalar quantum field theory after a quench from an equilibrium 
disordered initial state at the temperature T{ = T to a final state at Tf « 
jl7j . For unstable modes the instantaneous quench will be mimicked by a 
time dependent mass 

m 2 {t) =m 2 0{-t) -m 2 f O{t), 
which is tachyonic at t > 0. In the decomposition 

dk 1 

the corresponding mode functions Uk{t) obey 

d 2 



0(f ' t)= I (2% 2 PT- + a + (k)ul(t)e~^) , (32) 



dt 2 



+ k z + m\t) 



u k (t) = 



and initially (for t < 0) we have Uk(t) = e luJkt , Uk = \J m 2 + k 2 . For t > 
the solution is 

u k (t) = A k e akt + B k e- ak \ (33) 

with a k = J m 2 — k 2 (we will concentrate on unstable modes so that k 2 < 

m 2 ). Matching the t > and t < solutions and their first derivatives at 
t = 0, we can determine the A k and B k coefficients 

A k = Ul-i^), Bk = \(l + i ^). (34) 

The information about the domain size is encoded in the equal time 
correlation function (spatially averaged over the volume V) 

G(x, t) = i J dy<(f){x + y, t)<f>(y, t) > th . (35) 
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Indeed, if x is not greater than the domain size Lp, then + t)(j)(y, t) « 
(j) 2 (y,t) and the integral (ESI) should be near its maximal value G(0,t). On 
the contrary, if x ^> Lp, then the integral (ESD averages to zero. Substituting 
into ( 1351) and remembering that for our approximations 

< aa > th ~ 0, < a + a + > th & 0, 



we get 



1 



G(x,t) = 



dk 



V J 2u k 



u k (t)\ 2 < a(k)a + (k) > th e lk '*+ < a + (k)a(k) > th e 



-ik-x 



1 



dk 



V / 2ui 



u k (t)\ 2 < a(k)a + (k) > th + < a(-k)a + (-k) > th +6(0) 



(36) 

To understand the meaning of £(0), let us return one step backward and 
write 



KM I' 
2u k 



5(6) 



5(k - k') e-W*-*) ^0= dy = 

(2n) 6 y/2u k V2^fc' 



KMf 



\u k (t)\ 2 V 



2u k J (2tt) 3 2u k (2tt) 3 ' 
Besides, as was explained above, 

N 



< a + (k)a(k) > 



th z 



N = V/(2ir) 6 . 



Therefore, (ESI) takes the form 



G(£ >') = /(lFi K(t)|2cot H#J e * 

After the quench T « and hence coth (^) « 1. To study the growth of 
domains, one should subtract the contributions that were already present 
before the quench [47J 

r i 

G(x,t) = G(x,t) - G(f,0) = / t^t^- [\u k (t)\ 2 - 1] e ik ' x . 



(2tt) 3 2co k 



But 



and 



\u k (t)\ 



imply 
1 



l + cosh(2a fc *) + i (l 



tot 



a. 
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Therefore, the contribution of unstable modes in the domain growth is con- 
trolled by 

^ <> =\ I w? k i 1 + i) [cosh {2akt) - 11 e " = 



47r ^2 / ( &r ) sin ( &r )^ ( x + ^l) sinh2 dk > ( 37 ) 

o fc 

where r = \x\ and A: = \k\. We are interested in the large t asymptote of 
this function. Then 

sinh 2 (a k t) -> e 2c ^ 

and 

2 2 

G(x,t) ► — — 7T-J- I dk k— e 3 ^ \ 
167T z r z J oj k ct L k 
o 

where 

g(k) = 2a k t + In (At). 
Note that g(fc) has a maximum at ko « y^W an d 

m; 

Therefore 

2i 

p(/c) « g(ko) (k- k ) 2 

m f 

and for large t the function e 9 ^ has a very sharp peak at k = ko. Near this 
point (cjfctt 2 ,) -1 i s a slowly varying function. Therefore 

m,f 

6(3 , t) -» -A- e 2 ""< /, J = / dfc <*-*»>" sin (fcr) , 



7 

J o 



where we have used 

1 1 



-~ , fi'(fco) ~ 2m d + In (&or). 

^fc «fc mm} 



The remaining integral 



oo 



/ dke l( - k fco)2 sin (At) = sin (/c r) e *r, 7 = 



7 mj 

—00 
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As we can see gZ| HE] 



Therefore, the domain size grows with time, according to the Cahn-Allen 
scaling relation [48, 49j 

L D (t) = (38) 

V m f 

Remembering our estimate ( 1311 ) for the rolling-down time and taking a mean 
value of ( 1231 ) as an estimate for m!j: mj = |A s t> 2 , we get 

L ^ = -\/r^^ L4fm \A- ( 39 ) 

v V A s ir V v z 

Assuming Gaussian initial fluctuations of the fields and of their derivatives, 
it can be shown JIB] that the probability of the initial strength b of the 
axial-vector current to be large is exponentially suppressed. Therefore, our 
estimate ( 091 ) shows that typically DCC domains are quite small, unless b is 
large enough in some rare occasions. One concludes that "formation of an 
observable DCC is likely to be a rather natural but rare phenomenon" [50J. 

To close this section, let us indicate some review articles about the DCC 
phenomenon [40, H21 GHIl EI] , where an interested reader can find further 
discussions and references to the original literature, which is quite numerous. 

4 DCC at a photon collider 

Usually DCC formation is considered in the context of heavy ion or hadron- 
hadron collisions. We see no reason why gamma-gamma collisions have 
to be discriminated in this respect. The basic interaction for the photon 
is of course an electromagnetic interaction with charged particles, and to 
lowest order in a the photon appears as a point-like particle in its interac- 
tions. However, according to quantum field theory, the photon may fluc- 
tuate into a virtual charged fermion-antifermion pair. In this way strong 
interactions come into play through quark-antiquark fluctuations. While 
the high-virtuality part of such quark-antiquark fluctuations can be calcu- 
lated perturbatively, the low-virtuality part cannot. The latter is usually 
described phenomenologically by a sum over low-mass vector-meson states 
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- the Vector Meson Dominance (VMD) ansatz. Therefore, the effective 
photon state vector has the form [H21 53J 

|7 >= \/%\jbare > + ^c v \V > + ^c q \qq > +^q|/ + r >, (40) 

V q I 

where [52j 

4-^, £-2.2, |, 2 , 6 , |, 18 ,, ^,n,. 

The coefficients of the perturbative |gg > part depend on the scale \i at that 
the photon is probed. Namely [51^ 



2 _ a 21 

c„ e„ in « 5 

where |e g | = |, | is the quark charge and /co is an unphysical parameter 
separating the low- and high-virtuality parts of the quark-antiquark fluctu- 
ations. 

The last term in (JJQj) describes fluctuations into lepton pairs and is unin- 
teresting inasmuch as the hadronic final state is concerned. The coefficient 
of the first bare-photon term is given by 



V q I 

and is close to unity. 

Therefore, for some fraction of time the photon behaves like a hadron. 
This fraction is quite small, about 1/400 [54J, but for hadronic final states 
this smallness is overcompensated by the fact that in its hadron facet the 
photon experiences strong interactions. 

According to (110j) . one has six different types of possible interactions in 
the high-energy photon-photon collisions [52j: 

• both photons turn into hadrons (vector mesons) and the partons of 
these hadrons interact with each other. 

• one photon turns into a hadron and its partons interact with the quark- 
antiquark fluctuation of another photon. 

• both photons fluctuate perturbatively into quark-antiquark pairs and 
subsequently these fluctuations interact with each other. 
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• a bare photon interacts with the partons of the hadron which another 
photon was turned into. 



• a bare photon interacts with the quark-antiquark fluctuation of another 
photon. 

• bare photons interact directly in a hard process. 

In fact, in the total hadronic cross sections, the first two event classes dom- 
inate, the bulk of the contribution coming from the p°p° component of the 
first class |52j. Therefore, high-energy photon-photon collisions are very 
much similar to the hadron-hadron collisions, and if DCC can be formed 
in the latter case, it will be formed also in the former case. Of course, the 
photon-photon cross section is strongly reduced compared to the hadron- 
hadron cross sections (about 10 5 times). However, it is not improbable 
that a great deal of this smallness is overcome by somewhat more favorable 
conditions for the DCC formation in the gamma-gamma collisions than in 
the proton-(anti)proton collisions. The argument goes as follows. As has 
been mentioned above, for boost-invariant initial conditions, when the field 
depends only on the proper time r = \/t 2 — x 2 , the d'Alembertian equals 



The second term describes the decrease of energy in a covolume due to 
longitudinal expansion and brings an effective "friction" , which is necessary 
for quenching, into the equation of motion (101 |55j. The transverse (D = 
2) and spherical (D = 3) expansions can be modeled analogously if one 



Therefore, the larger is D the more efficient is the quenching and the spher- 
ical expansion seems to be the most favorable for pion zero mode amplifica- 
tion [55IEB]. This simple observation is confirmed by a more detailed study 
|56j . However, to organize an isotropically expanding fireball is not a trivial 
task even in head-on hadron-hadron collisions. Constituent quarks inside 
hadrons become "black" at high energies, and for the projectile remnants 
not to spoil the isotropic expansion, one may wait for a rare event when 




assumes that the field depends only on r 




□ = 



dr 2 



D d_ 

t dr 
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these black disks inside the projectiles are aligned [S7j. The probability 
that all six constituents are aligned in colliding protons during a head-on 

collision is pi ~ f , while the analogous probability for the four con- 
stituents of p°p° collisions is P2 ~ \ ■ Taking r q « \r p « |r p , we get 

Vilvi ~ 10 3 - Therefore, photon-photon collisions seem to be more favorable 
in this respect. 

Even if the "right" fireball is prepared, the odds of the large DCC do- 
main formation are usually small. In [HE] this probability was found to be 
about 1CT 3 . Remember that our preceding considerations indicate that one 
needs large initial strength of the axial- vector ,S£7(2)-current. In gamma- 
gamma collisions this initial strength is expected to be enhanced due to 
chiral anomaly effects, analogous to what was considered in [58j for heavy- 
ion collisions. The effects of chiral anomaly can be incorporated in the linear 
sigma model by adding the following interaction Lagrangian 

£ - a 7T e F [W F aT - — ir°E ■ H 

t-'anom. — . f c fii'aT- L 1 — „ n J - J 11 • 
Jit ^ Jit 

Under SUa(2) transformations (JToT). we have 

5£>anom. = ~ E • H O 63. 
KJjt 

Therefore, according to the Gell-Mann-Levi equation (we neglect explicit 
symmetry breaking): 

and the axial-vector current is no longer conserved even in the limit of 
zero quark masses. As a result, electromagnetic fields can lead to the en- 
hancement of the axial-vector current strength in the 7r°-direction. The 
corresponding induced strength equals 

h = -4- / e ( t ) ■ H ( r ) < t ) rdr - ( 4l ) 

As was shown in the expected effects are small in relativistic heavy- 
ion collisions, but nevertheless this initial small "kick" can have substantial 
effect on the DCC formation. Unfortunately we can not use (14 ip to estimate 
how big is the kick in gamma-gamma collisions - the application of the sigma 
model makes sense only after some time after the collision, while the effects 
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of the chiral anomaly on the 63 magnitude are confined to the first instants 
of the collision. 

It was suggested f21 ED] that in hadron-hadron collisions DCC could 
be formed through the "Baked Alaska" scenario. However, as we have 
mentioned above, the considerable part of the 77 — > hadrons cross section 
is due to the p°p° mechanism. Therefore, the Baked Alaska model should 
work for gamma-gamma collisions too. So let us take a closer look at it. 

Normally Baked Alaska is a delightful dessert where ice cream is covered 
by meringue and then baked very quickly in a hot oven without melting the 
ice cream (you can find the recipe in the appendix. Try it and enjoy). In 
physics context, however, this term firstly appeared as denoting a model for 
nucleation of the B phase of superfluid 3 He inside the supercooled A phase 
|nH E21 E3j- The surface tension at the boundary between the A and B 
phases is anomalously large. Therefore usual small bubbles of the B phase, 
created inside the A phase by thermal fluctuations, are energetically not 
profitable. Hence, they shrink and vanish. Only for a very large bubble the 
volume energy gain overcomes the surface energy and the bubble begins to 
grow. However, it is virtually impossible to create such a gigantic critical 
bubble by thermal fluctuations. The experiment, nevertheless, discovered a 
high enough nucleation rate. To explain the puzzle, Leggett suggested that 
the nucleation was assisted by cosmic rays |HT]. Secondary electrons from 
the passage of a cosmic-ray muon through the liquid create hot spots in 3 He 
by depositing several hundred eV energy in small volumes. Inside such a 
"fireball" the Cooper pairs of the S H e atoms are broken and, therefore, the 
normal, Fermi-liquid phase of the s He is restored. Yet, the fireball expands 
quickly and becomes a "Baked Alaska": a cold core surrounded by a hot, 
thin shell of normal fluid. There is some probability that after the core is 
cooled below the superfluidity phase transition temperature, it finds itself 
in the B phase. This B phase core can expand to larger than the critical 
radius because for some time it is shielded from the A phase bulk by the 
expanding hot shell, thereby eliminating surface energy price of the A — B 
boundary layer. When the shielding shell finally disappears, the B phase 
bubble is larger than the critical one and, therefore, expands further until 
it fills the whole vessel. 

Let us return to Baked Alaskas produced by high-energy gamma-gamma 
collisions. Suppose DCC is formed inside the Baked Alaska core with the 
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misalignment angle 6. That is inside the DCC region one has 

< CT >DCC= fn COS 6, < 7T >DCC= fn sin & ™: 

n being a unit vector in isospin space. Outside the fireball one has the 
normal vacuum: 

< a >= f n , < 7f >= 0. 

Finally, when the shielding shell of hot hadronic matter disappears, the 
DCC relaxes to this outside normal vacuum by emitting coherent low energy 
pions. Hadronization of the shell also produces mainly pions and, therefore, 
generates a background to the DCC signal. Simple considerations allow one 
to estimate the numbers of the DCC and background pions jSO] . Energy 
density in the DCC region is higher than in the normal vacuum because of 
the symmetry breaking term Vsb = —Ha. The difference is 

Ae v = -H<a > DCC +H < a >= Hf^l - cos 0) = 2/> 2 r sin 2 ^. 

Therefore, the total volume energy available for pion radiation from the 
DCC decay is 

Ev = y^/X sm 2 ' 
where R is the fireball radius at the moment of hadronization. However, 
pions radiated from the DCC are nonrelativistic (in the DCC rest frame). 
Therefore, the expected average number of such pions is (we have assumed 
that < sin 2 f >= \) 

NyM^^tffim,. (42) 
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One can assume [60] that at the moment of hadronization the shell consists 
of one densely packed layer of pions, each having the radius r n « \m~ l « 
0.7 fm. Therefore, the number of background pions from the fireball shell 
is 

N b « — . (43) 

For a large DCC bubble of the radius R « I0r n w 7 fm, the above given 
estimates imply 

N v « — 125 (—} w 250, N b w 400. 

3 V m 7r/ 

There will also be coherent pions associated with the surface energy of the 
interface between the DCC and outside vacuum. The energy density in the 
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interface is dominated by the contribution due to gradients of the fields. If 
one assumes that the interface thickness is the same as for the hadronized 
shell, that is d ~ 2r w « wi" 1 , then 



e s « \ [(Att) 2 + (Aa) 2 ] « ^ [sin 2 # + (1 - costf) 2 ] = 2§ sin 2 6 



2 2d 2 v 7 J d 2 

The corresponding total average energy is thus 

E s « 4tt R 2 d < e s >« 47T i? 2 / 2 m^. 

Pions originated from the surface layer of thickness <i will have characteristic 
momenta ~ 1/c? « and energy i?^ = + « ^/2m n . Therefore 
the expected average number of surface pions is 

N s « — « 2V2ttR 2 f 2 « 50\/27r f — ^ « 105. 

As we can see, the DCC signal from such a large single domain is quite 
prominent. This becomes especially evident if one realizes that there is a 
large probability, ~ 10%, that almost all of these 250 nonrelativistic signal 
pions are charged ones, with only a few neutral pion admixture. The prob- 
ability that such a huge isospin-violating fluctuation happens in the back- 
ground pions is, of course, completely negligible. This striking feature of 
the DCC signal follows from the following simple geometrical argument p3! 
(the inverse-square-root distribution, discussed below, was independently 
rediscovered many times by different authors. See jlOj for relevant refer- 
ences). Let the pion field in a single DCC domain be aligned along a fixed 
isospin direction n = (sin 9 cos 0, sin 9 sin 0, cos 9) . Classically the radiation 
is proportional to the square of the field strength. Therefore the fraction of 
neutral pions /, emitted during relaxation of such a DCC domain, equals 

I |2 

/ = J 7 " 31 = cos 2 9. (44) 

E kl 2 

a=l 

The probability for / to be in the interval (/, / + df) is given by 

P(f)df = [P(cos 9) + P{- cos 9)} d cos 9. (45) 

Any orientation of the unit vector n is equally valid. Therefore the proba- 
bility P(cos9)dcos6 for finding cos 9 in the interval (cos 9, cos 9 + dcos9) 
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equals 

cos 0+d cos 8 2tt 
If f 1 

— / dcos9 / dd> = - d cos 9 . 
4tt J J Y 2 

cos6> 

This implies P(cos6) = |. Therefore, from ( 1451) the probability density for 
the neutral fraction / is 

PU) = 

while from (JUJ) c/f = 2 cos d cos 6 = 2v7 dcosO, and we finally obtain 

- (46) 

This inverse-square-root distribution is drastically different from what is 
expected for noncoherent pion production: the binomial-distribution which 
for large pion multiplicities TV turns into a narrow Gaussian centered at 
f = h 



For example, the probability that / does not exceeds 0.01 according to (J3Ej) 
equals 



o.oi 



p(/<o.oi) = j ^= = Vom = io% !!, 



while the binomial-distribution predicts 

o.oi 

F(/<0.01)= J P nc (/)4fw0.01^e-T<l, foriV»l. 
o 

However, as we have seen above, it is not easy to produce large DCC 
domain. For many small DCC domains, with random vacuum orientations, 
the effect of the inverse-square-root distribution will be washed out by av- 
eraging over the orientations (it was, however, argued in [65j that the later 
case of small DCC domains may lead to enhanced baryon-antibaryon pro- 
duction within the framework of the Skyrmion picture of the nucleon). The 
intuitive reason why it is difficult to grow up a large DCC bubble is the 
following. Up to now the only mechanism discussed by us for the DCC 
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formation, was the spinodal instability which operates when the fields are 
rolling-down from the top of the Mexican hat potential towards the valley. 
But this rolling-down time is small unless the initial strength of the axial- 
vector current is large (see Eq. I3T1) . Fortunately, there exists a parametric 
resonance mechanism (HE) which can further assist the DCC formation, af- 
ter the spinodal instability is over. To illustrate the physical idea, let us 
again return to the Blaizot-Krzywicki model jlBj with boost-invariant initial 
conditions. The equation of motion for the pion field is 

□tt = -A s (a 2 + 7f 2 - v 2 ) 7T (47) 

with 

D ^ d 2 Id 
dr 2 t dr 

At later times the fields are near the true vacuum. So we assume 

v ~ A, 7? = 7r(r)n. 

The above pionic field describes a (small) disoriented chiral condensate 
aligned along a fixed unit vector ft in isospace - the result of preceding 
spinodal instability. If we neglect nonlinear terms, we get in the zeroth 
approximation (note that A s (/ 2 — v 2 ) = m 2 ) 



7r(°> (r) + ^-^ + ml 7r<°> (r) = 0. (48) 



*(°) 
r 

This equation is equivalent to the Bessel equation and its general solution is 
a linear combination of the Bessel functions Jo(m n r) and Y^m^r). There- 
fore in the large r limit one expects damped oscillations 

7r(°Vr) = —= cos (m^r + (/?), 

with A and <p as some constants. Now consider a fluctuation (in the re- 
direction) 7r^(r) around the zero-mode 7P°)(t): 7t(t) = 7P°)(t) +7P 1 )(t). 
Keeping only the terms linear in tt^\t) we get from (JUJ) and 



TfW(r) + + ml TrW(r) = -3X s tt^ 2 (t) 7t«(t), 

r 



or (we have neglected unimportant 3X ° A tt^(t) term) 



7f (1) (r) + ^^ + a; 2 1 + 1 cos (lot + 2<p) tt (1) (t) = 0, 



<7 
r 
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where loq = m n , u = 2171^ and q = 2 S 2 . Therefore, one expects a paramet- 
ric resonance because to = 2coq. 

For more general initial conditions, parametric instabilities are expected 
for the low momentum pion modes The energy of the cr-field oscil- 

lations around a = can also be pumped into pionic modes through 

the parametric resonance |66[ 158] . In this case, naively only modes with 
r~2 

k ~ yj -f — ml « 270 MeV can be amplified, because oscillation frequency 
in the cr-direction is m a = 600 Me^. However, the energy can be redis- 
tributed in the long wavelength modes due to nonlinearity. 

An interesting example of the parametric instability is Faraday waves 
|69| - surface waves parametrically excited in a vertically vibrating con- 
tainer of fluid when the vibration amplitude exceeds a certain threshold. 
The resulting standing waves on the fluid surface can form funny intri- 
cate patterns jZQj. Even a more closer analog is given by the induction 
phenomenon in quartic Fermi-Pasta-Ulam (FPU) chains |71j : the energy, 
initially supplied to a single harmonic mode, remains in this mode over a 
certain period, called the induction time, when it is abruptly transferred 
to other harmonic modes. The original explanation [7T] involves nonlinear 
parametric instabilities similar to the one described above. It should be 
mentioned, however, that our above arguments, in favor of the exponential 
growth of fluctuations due to the parametric resonance, are heuristic. A 
naive perturbation theory, implicit in these arguments, is not adequate for 
such nonlinear problems. In the case of the FPU chains, it was argued that 
a more correct treatment was provided by shifted-frequency perturbation 
theory [72J or by Krylov-Bogoliubov-Mitropolsky averaging technique |73j. 
We are not aware of similar studies in the context of DCC dynamics, but 
the reality of parametric instabilities is indirectly confirmed by numerical 
studies [34, 74J: the observed amplification of long- wavelength pionic modes 
last much longer than expected solely from the spinodal instabilities. Be- 
sides, the amplification of pionic modes with k « 270 MeV was clearly 
demonstrated. 

5 Quantum state of the disoriented chiral condensate 

The eventual decay of the DCC is a quantum process, because one registers 
pions and not the classical field. Therefore, the natural question about the 
DCC quantum state \ij >dcc arises. The usual way of quantizing some 
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classical field configuration is to use coherent states which are eigenstates 
of the annihilation operator JZHJ 

a \a >= a \a > . (49) 

Decomposing 

\a >= 2_^c n \n >, \n >= ' |0 > 



n=0 

and using a\n >= y^n\n — 1 >, we get from the recurrent relation 
\/nc n = ac n _i. Therefore 

a n co 

71 j 

and 

\a >= Co exp (aa + ) |0 > . 

However, 

< a\a >= < a|n >< n\a >= |c n | 2 = |co| 2 exp (a* a), 

n n 

and, therefore, the normalization condition < a\a >= 1 determines cq up 
to a phase. Finally 

— + aa + I |0 > . (50) 

The generalization of this construction to the DCC classical field configu- 
ration f(x) is [7Bj (the isospin indices are suppressed) 

\r) >dcc= exp (-1 J dkf\k)f(k) + J dkf(k)a + (k)^j |0 >, (51) 

where f(k) is the Fourier transform of the DCC classical field. 

However, it was argued [HI EE] that ^ ne quantum state of the disoriented 
chiral condensate is expected to be a squeezed state , if the parametric 
amplification mechanism, discussed above, indeed plays a crucial role in the 
DCC formation. To explain why, let us consider a one-dimensional, unit 
mass, quantum parametric oscillator with the Hamiltonian 
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where the p and x operators are time-independent in the Schrodinger picture 
and obey the canonical commutation relation (K = 1): 

[x,p] = i. 

Quantum state vector \ip > of this oscillator is determined by the Schrodin- 
ger equation 

ij t \^>=H(t)\^> . (52) 

Lewis and Riesenfeld gave [80j a general method of solving the Schrodinger 
equation by using explicitly time-dependent invariants which are solutions 
of the quantum Lieuville-Neumann equation 

§ + i[H, I] = 0. (53) 

It turns out that the eigenstates of such a Hermitian invariant I(t) are just 
the desired solutions of the Schrodinger equation up to some time-dependent 
phase factor. Let us demonstrate this remarkable fact jSDj • The eigenvalues 
of the Hermitian operator I(t) are real. Therefore, from ( 1531) one easily gets 

r)T 

i < A'| j=- |A >= (A - A') < X'\H(t)\X >, (54) 

where |A > is an eigenvector of the operator I(t) with the eigenvalue A: 

I(t) |A >= A |A > . (55) 

In particular 

< A| g |A >= 0. (56) 

By differentiating (1551) with respect to time and taking the scalar product 
with | A' >, we get 

< A'| — |A >= (A - A') < A'| - |A > + 5\\i — . (57) 
For A = A', we get from ( 1371) and 



<9A ^ 

ai = °- 

That is the eigenvalues of the operator I(t) are time-independent (as it 
should be for the invariant operator). From ( 1371) and ( 1541) one gets 

(A - A') < A' | £ |A >= (A - A') < \'\H(t)\\ > 
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and, therefore 



< A'| £ |A >=< X'\H(t)\X >, if X' ^ X. (58) 



If (EBD is also satisfied for the diagonal matrix elements (A = A'), then we 
can immediately deduce that |A > is a solution of the Schrodinger equation. 
But this may not be the case for our particular choice of eigenvectors. Nev- 
ertheless, in this case we can still adjust the phases of the eigenvectors in 
such a way that the new eigenstates 

|A >'= e i6x[t) \X >, A (O) = O, 

insure the validity of (EED for all A, A'. All what is needed is to choose the 
time-dependent phases 6\(t) in such a way that one has 

< A| e"** ( i^- ) e** |A >=< X\H(t)\X >, 



or 



Therefore 



dt 



f-<A|i|-A W |A>. 

t 

0x(t)= I <X\i^-H(r)\X> dr. 



To summarize, we can take any set |A > of the eigenstates of the invariant 
operator I(t) and express the general solution of the Schrodinger equation 
as a linear combination 



t 



|^>=^C A exp|z J <X\z^--H(r)\X> dr^\X>, (59) 

with some time-independent coemcients C\. 

To find the invariant /(t), let us firstly construct its classical counterpart 
I{t) by the simple and transparent method of Eliezer and Gray [EI]- The 
equation of motion 

x + u 2 (t) x = 

can be viewed as the ^-projection of a two-dimensional auxiliary motion 
governed by the equation 

f+u 2 (t)r = 0, (60) 
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where r ' = xi + y j . In the polar coordinates 



x = pcosp, y = psmp, 



and 



r= (p- pip ) e p + {pp + 2pp)e (p , 
where the unit basic vectors are 

e p = cos p i + sin pj, = — sin p i + cos ^ j. 

Therefore, ( 1601) in the polar coordinates takes the form 

1 d(p 2 p) 



p — p<p +u (t)p = 0, p{p J r2p(p = - 



p dt 



0. 



The second equation implies 



p ip = L = const. 



This is nothing but the conservation of the angular momentum for the 
auxiliary planar motion, when the first equation can be rewritten as 

T 2 

p + u (t)p = 



P 



3 ' 



Let us now remark that (for unit mass) 



p = x = p cos p — pp sin p 



px — L sin p 
P 



Therefore, 



and 



L 

L sin p = px — pp, L cos p = —x, 

P 



L 2 

-^-x 2 + (pp — px) 2 = L 2 = const. 
P 2 

The Ermakov-Lewis invariant JH21 [83J corresponds to the particular case 
when the angular momentum L has a unit value 

,2 



/(*) = \ 



X 



^ + {PP~ px)' 



p* 



It is straightforward to check that its quantum counterpart 



i(t) 



i 

2 



x 



.P' 



+ {pp - p x Y 



(61) 
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obeys the quantum Lieuville-Neumann equation ( 1531 ) if the auxiliary func- 
tion p{t) satisfies the Ermakov-Milne-Pinney equation jB3j 



p + uj 2 {t)p = — 



P" 



(62) 



To find eigenvectors of the operator I(t), let us note that 

/(*) = b + (t)b(t) + \ , 



(63) 



where we have introduced time-dependent "creation" and "annihilation" 
operators 



b(t) = 



i 



V2 



x 



p 



+ i(pp - px) 



b+it) = 



1 



V2 



x 



P 



i{pp — px) 



(64) 



It can be immediately checked that one indeed has the canonical commu- 
tation relation 

[6(t),fc + (t)] = l. (65) 

Equations ( 1531) and ( 1551) indicate that the eigenvectors of I(t) are 6-number 

states 

(i+(t))> 



n; 6 >= 



0:6 >. 



where the 6-vacuum state is defined by the condition 

bit) |0;6 >= 0. 



(66) 



(67) 



Using these eigenvectors, one can construct the solution of the Schrodinger 
equation, as described above. In the DCC case, however, it is preferable 
to express the quantum state vector in terms of the pion creation and 
annihilation operators a + , a. So we need a relation between two sets of 
the creation-annihilation operators a + , a and 6 + , b (the Bogoliubov trans- 
formation). The pionic modes correspond to the late time asymptotes 
— ► cj(oo) = cjo- The solution of the Ermakov-Milne-Pinney equation 
for u(t) = luq = const is 

1 

P = 



UJ 

Therefore, the "pionic" creation-annihilation operators have the standard 
form 

1 a a 1 

a = [loq x + ip], a + = [cjq x — ip] . 



(68) 
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The comparison of ( 1681) and 
tion 

b(t) = a(t) a + p*(i) a+, b + (t) = (3(t) a + a*{t) a 
where 



gives the desired Bogoliubov transforma- 

(69) 



a(t) = 

(3{t) = 
It can be checked that 



1 



1 



1 



At) 
i 

W) 



+ u p(t) -ip(t) 



u p(t) +ip(t) 



(70) 



|a(t)| 2 -|/?(t)| 2 = l. 
Therefore, up to an irrelevant common phase, we can take 

a{t) = coshr(t), (3{t) = e i5 ® sinhr(t). 

The Bogoliubov transformation 
mation 

b = S{z)aS + {z), b + = S{z)a + S + {z 
where S(z) is the so-called squeezing operator [79j 

1 



(71) 

can be viewed as an unitary transfer- 

(72) 



S(z) = exp 



( 



za a 



y* aa) 



z = r e 



i(6+ir) 



(73) 



1 



Indeed, using the Campbell-Hausdorf formula 

eHe-z = A + [B, A] + ±[fl, [5, i]] + -[5, [5, [5, A]]] + 

and summing up the resulting infinite series, one can check that ( 1721 ) and 
(ESI) are equivalent, if a and (3 are given by (EH). 

Therefore, the DCC quantum state is expected to have the form 

\V >DCC= s(z)\v Q >, 

where > is the DCC initial state before the onset of the parametric 
amplification. If |\I/o > is a coherent state, then the resulting \r] >dcc state 
will be the one called the squeezed state [721 • The actual parameters of this 
state (for example z) are determined by the initial conditions and are hard 
(if not impossible) to estimate from the theory alone. 

Let us now recall the isospin, for a moment. The isospin generators are 



1=1 7?(x) x tt{x) dx. 
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The classical pion field 7r of the DCC points in some fixed direction, n in 
the isospace. If its time derivative 7? also points in the same direction then 
1 = and we will have an isosinglet state. One can expect such a situation 
in p — p collisions (which dominates in the 77 — ► hadrons cross section), 
because the "vacuum cleaning" effect, which precedes the DCC formation, 
in this case is mainly due to two colliding isospin blind gluon walls. It is 
easy to construct an isosinglet squeezed state by just exponentiating the 
apparently isoscalar operator 

3 

— af a\ = 2a+a~t — a^a^ , 

i=l 

where a± = ^iflt ± i^t) are the charged-pion creation operators. The 
resulting squeezed state is j7S] 



I* >= N exp || (2a%at - ajaj) } |0 > . 

Well, this expression does not seem, at first glance, to correspond to the 
canonical form ( 1731) of the squeezing operator. In fact, it indeed gives a 
squeezed state. This is clear from the following normal-ordered form of the 
squeezing operator [79j 



Sire 49 ) = N exp (| a+a+) ~ H ~ \ a+ Ti a T ex P ("f aa ) ' ( 74 ) 
where a = e i6> tanh r and 



iv = ^^ = (i-H 2 ) 1/4 . 

V cosh r 

What is the probability P(m; n) that > decays by producing in total 
2n pions and among them 2m neutral pions and equal numbers of positively 
and negatively charged pions? According to the standard rules of quantum 
mechanics 

P(m; n) = \ < m; n|\I/ > | 2 , 
where the normalized state \m;n > is defined through 



\m; n >= — = — — (atat) [n ~ m) (at) 2m |0 > 

However, 



\^>=NJ2 {2a + + a + _ - a+a+) k |0 >= 

k={) 
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3C 



= N E ^~jr~ E c i 2 " (J - k)\ \l; k > 



fc=0 



fc! 



Z=0 



Therefore, 
P(m: n) 



N ^-j- C™ 2 n - m 



m 



^/(2m)\ 



in — m)\ 



= N 2 \a 



2n 



(2m)! 
(ml 2 m ) 2 



One can prove by induction in n that 

(2m)! (2n+l)! 



E 

771=0 



m 



n 



\2 n ) 2 



This enables one to express P(m; n) as a product of two probabilities: 

P(m; n) = Pi(n)P2(m; n), 

where Pi(n) is the probability that one will find the total number of 2n 
pions after the state |^ > decays 

P l( n) = ^H-^ +1)! 



(n! 2-) 2 



Note that 



E 

71=0 



(2n + l)< |2n 



! 2") 2 



3-5 
2~ 4 



3-5-7 
2-4-6 



+ = (l-\a\ z ) 



2\-3/2 



Therefore, the normalization coefficient N should be N = (1 — |a| 2 ) 3 / 4 , and 
this is exactly what is expected from ( 1711) for the product of three properly 
normalized single-mode squeezed states of Cartesian pions. 

More interesting for us is the second factor, P 2 (m;n), the probability 
that one finds 2m neutral pions in such a 2n-pion final state JHDl E5] : 

[n\2 n ) 2 (2m)! 



Pofm; n) = 



(2n + 1)! (m!2 m ) 2 ' 

In fact, P 2 (m; n) is a particular case of the Polya distribution 
n are both large, one can use the Stirling formula 

n! ~ ( — 

V e/ 

to get 

1 1 

P(m; n) ~ — 



If m and 



irn 



2n ^/m~J 



n 
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Therefore the same inverse-square-root distribution ( 1461) is recovered in the 
continuum limit. 

A few more words about coherent and squeezed states, in order to make 
them more familiar. If [A, B] is a c-number, then [87j 



Mb = e A e B e -\ [A,B] 



Using this theorem, we get 



exp [a a —a a\ = exp 



1 * 
-a a 



exp (aa + ) exp (—a* a). 



But exp(— a*a)|0 >= |0 >. Therefore, the coherent state \a > can be 
generated by the unitary displacement operator: \a >= D(a) |0 >, where 



D (a) = exp [a a + — a* a}. 



(75) 



Now let us consider the ground state for a harmonic oscillator (we have 
abandoned the unit mass restriction but will still keep h = 1): 



^o(x) =< x\0 >= [27R7o] 1/4 exp 



X 

2^ 



and calculate the effect of the displacement operator on it. For the harmonic 
oscillator 



1 



a = 



V2 



[mcu x + i p] , a 



1 



V2 



mu 



mcu x — ip]. 



Hence, 
where 



a a —a a = ipo x — ixq p, 



xo = \ Re a, p$ = \/2mu Im a 

V mid 

In the coordinate representation p = —i-§-] therefore , 

d 



ty cs (x) =< x\a >= exp 



IPQX — Xq 



dx 



= [27rcro] exp 



x - Xq\ .xopo 
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As we see, for harmonic oscillator, the coherent state is a Gaussian which 
is displaced from the origin by xq. It has the ground state width gq and 
a phase linearly dependent on the position x. Schrodinger discovered |BBj 
such a state as early as 1926 while seeking "unspreading wave packets" . 

What about the squeezed vacuum state ^ s o = S(z) ^o? Note that (for 
simplicity we will take z = r to be real) 



a + a + — aa = — i(xp + px). 



Then 



-it 



xp + px, X 



-IT 



■rx, 



2 ' 2 

and the Campbell-Hausdorf formula will give 



xp + px, p] = rp, 



S(r)xS + (r) = e r x, S(r)pS + (r] 



T ~ 

e p. 



Therefore 



9r ~9 9 — 9r * 9 



1 



<9 2 



2m ' 2 2m d{e- r x) 2 ' 2 

and comparing the equations which ^ s q and are satisfying 

we conclude that 

V s0 (x) = ^o{e- r x) = [2ttcjIY 1/A exp 



+ — — (e r x) L 



x 



2e r a 



In a more general case, the squeezed state ^ ss = D(a)S(r) is a Gaussian 

2 



^ ss {x) = [2tt<7§] 1/4 exp 



X — Xq 

2a 



+ ipox 



x po 



where 



a = e ctq. 



Therefore, the squeezed state differs from the coherent state only by squeez- 
ing in (if e r < 1), or squeezing out (if e r > 1) the width of the ground state 
Gaussian. The momentum-space wave function is squeezed oppositely. 

For further discussion see O EH| • Let us note that squeezed states were 
also discovered long ago (in 1927) by Kennard [87J. Both the Schrodinger 
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and Kennard works had little impact, and both coherent and squeezed 
states, which are now cornerstones of quantum optics, were rediscovered 
after decades. As Nieto remarks, "To be popular in physics you have to 
either be good or lucky. Sometimes it is better to be lucky. But if you are 
going to be good, perhaps you should not be too good" . 

As the last, but not the least, question of this section, let us ask whether 
a quantum state of the disordered chiral condensate \r) >dcc m &y be pro- 
duced without any intermediate phase transitions altogether, through the 
quantum reaction 77 — > \r) >dcc- The functional integral methods and 
stationary phase approximation (semiclassical approximation) are natural 
tools to study the scattering amplitudes between initial wave packet states 
and certain final coherent states fHH E] • In this article we are not really 
interested in actual calculations of this type. Our aim at the beginning is 
more humble - to provide some arguments that such a quantum transition 
is indeed possible and interesting. So we will consider an oversimplified toy 
model with the (second quantized) Hamiltonian 

H = to a + a + 2cu b + b + g (ba + a + + 6 + aa), 

with 2cj = m n . Here the 6-mode mimics neutral pions, a-mode - photons, 
and the interaction term with g = ^y 2 - imitates the 7r°77 interaction due to 
axial anomaly. Such Hamiltonians are used in quantum optics to study the 
second-harmonic generation [92, 93j. We further assume that all available 
initial energy, yfs, is accumulated in the a-mode. That is, the assumed 
initial a-mode occupancy is 

N _ y/s _ V£ 

It can be easily checked that [Ho, Hint] = 0, where 

Ho = cu a + a + 2u b + b, H{ n t = g (ba + a + + b + aa). 

This means that both //q and Hint are constants of motion. We can simply 
forget about Ho, because it just gives an irrelevant common phase factor 
exp (iN a u)i) in the evolution operator. Hence the nontrivial part of the 
initial-state evolution is given by the relation 

\^{t) >= e^ mtt \N a ;0> . 

However, N a ^> 1; therefore, as far as the 6-mode initial development is 
concerned, we can replace the a and a + operators in H int by the c-numbers 
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(76) 



As we see, a coherent state of 6-quanta is being formed. The mean number 
of quanta in this state grows with time (until the approximation considered 
breaks down) as follows: 



To estimate the terminal time, let us note that the variance of the number 
of quanta in the coherent state ( T7BI ) also equals |/3| 2 (see, for example, [89j ) . 
Therefore, the development time for the coherent state ( I7BI) can be estimated 
from the energy-time uncertainty relation 



Using this estimation, we finally get from ( 1771 ) the relation 



Therefore, for example, AT& ~ 100 implies yj~s ~ 2 TeV. 

Having in mind a very crude and heuristic nature of our arguments, we 
admit that we may easily be wrong by an order of magnitude. Nevertheless, 
the main indications of the above exercise that the axial anomaly can lead 
to the generation of a pionic coherent state in gamma-gamma collisions, and 
that the efficient generation requires not fantastically high center-of-mass 
energies certainly seems interesting and deserves further study. 

6 Concluding remarks 

"Have no fear of perfection - you'll never reach it" - Salvador Dali once 
remarked. At the end of our enterprise we reluctantly realized how true 
the second half of this quotation is. Therefore, we abandon an unrealistic 
dream to produce the perfect review of the disoriented chiral condensate 
and try to finish here. It is important to finish at right time, is not it? One 
American general began his talk with the following sentence: "My duty is 
to speak, and your duty is to listen. And I hope to end my duty before 
you end yours" . We hope that the reader has not yet end his duty, because 
there is one topic which should be touched a little before we finish. 



N b = \P\ 2 = (gN a t) 2 . 



(77) 
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The disoriented chiral condensate is a very attractive idea, and it has 
some solid theoretical support behind it, as we tried to demonstrate above. 
But are there any experimental indications in favor of its real existence? 
In fact there are some exotic cosmic ray events, called Centauros, where 
one may suspect the DCC formation. Centauros were discovered in high- 
mountain emulsion chamber cosmic ray experiments |94|. Typically, the 
detectors used in such experiments consist of the upper and lower chambers 
separated by the carbon target. Each chamber is a sandwich of the lead 
absorber and the sensitive layers. The normal cosmic ray event is usually 
generated by the primary interaction at about 500 — 1000 m altitude above 
the detector apparatus. About one third of the products of the primary 
interaction are neutral pions. Each neutral pion decays into two 7 quanta. 
Therefore, roughly one 7 quantum is expected per a charged particle in the 
primary interaction. When the interaction products reach the upper cham- 
ber the numbers of electronic and photonic secondaries are much increased 
through the electromagnetic shower formation. Therefore, the upper cham- 
ber usually detects several times more particles than the lower chamber, 
because the electromagnetic component is strongly suppressed by the car- 
bon layer, leaving mainly the hadronic component to be detected by the 
lower chamber. A big surprise was the discovery of events with the con- 
trary situation. Such events were named "Centauros" because it was not 
possible to guess their lower parts from the upper ones. 

The first Centauro was observed in 1972 at the Chacaltaya high mountain 
laboratory It was initiated by the primary interaction at a relatively 
low altitude, at only (50±15) m above the detector. Therefore, the event was 
very clean, that is was almost not distorted by electromagnetic and nuclear 
cascades in the atmospheric layer above the chambers. After correcting 
for the hadron detection efficiency and for the influence of the secondary 
atmospheric interactions, the event can be interpreted as the production 
of only one electromagnetic (e/7) particle and 74 hadrons with the total 
interaction energy ~ 330 TeV. 

Afterwards some more Centauros were found. Namely [94|, the Chacal- 
taya experiment observed 8 unequivocal Centauros, and two experiments at 
Pamir found 3 and 2 more Centauros. But no clean Centauros were found 
in Kanbala and Fuji experiments - the puzzle which still remains a mystery 
jni|. However, if the definition of Centauro is somewhat relaxed and all 
hadron-rich species are considered, then such Centauro-like anomalies con- 
stitute about 20% of events with the total visible energy > 100 TeV |94j, 
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hence they are by no means rare phenomena at such high energies. 

Of course the DCC formation is a candidate explanation of Centauro 
events. However such explanation is not without difficulties JS1]. For ex- 
ample, the large transverse momenta observed in the Centauro events is 
difficult to explain in the DCC scenario. It is not also evident that the 
Centauro hadrons are pions. If they are mostly baryons instead, then an 
alternative explanation may be provided by strangelets [95j . 

Let us note, however, that if one takes the DCC explanation of Cen- 
tauros seriously, some predictions immediately follow. First of all, it may 
happen that the DCC domain is produced in the cosmic-ray interactions 
with significant transverse velocity. In this case the coherent pions from 
the DCC decay will constitute "coreless jet" in the laboratory frame, with 
pions in the jet having small (< 100 MeV) relative transverse momenta [40J. 
Interestingly, such hadron-rich events, called Chirons, were really observed 
in both Chacaltaya and Pamir experiments [94|. 

If the DCC is aligned along the 7r°-direction in isospace, then a partic- 
ular anti-Centauro event is expected with neutral pion fraction / close to 
unity. For example, one has not very small probability that the neutral pion 
fraction is in the interval 0.99 < / < 1: 



No such anti-Centauro events were observed in the Chacaltaya and Pamir 
experiments. However, some anti-Centauros were reported in the Japanese- 
American JACEE experiment, in which the emulsion chambers were flown 
near the top of the atmosphere by balloons [94j . By the way, this experiment 
has not seen any Centauro events - another mystery puzzle of this cosmic 
ray Centauro business. 

Of course Centauro-like events were searched in accelerator experiments 
|94| . The first searches have been performed by UA1 and UA5 experiments 
at CERN even before the DCC idea was suggested. Both experiments found 
no Centauro candidates in the central rapidity region. 

The estimated average energy of Cosmic-ray Centauro events is about 
1740 TeV jp4j . If Centauros are formed in nucleon-nucleon collisions, this en- 
ergy threshold translates into y/s ~ 1.8 TeV in the cm. frame - roughly the 
Tevatron energy. This observation maybe explains the failure of UA1 and 
UA5 experiments, where the maximal available energy was yfs ~ 0.9 TeV, 



1.0 




0.99 
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and makes Fermilab experiments more attractive in this respect. However, 
one has to bear in mind that there is a crucial kinematic difference be- 
tween cosmic-ray and collider experiments [94] : the cosmic-ray experiments 
generally detect particles from the fragmentation rapidity region whereas 
the collider experiments are mainly focused on the central rapidity region. 
Therefore, the fact of observation of cosmic-ray Centauros does not auto- 
matically guarantee that these beasts can be found in Tevatron experiments. 

A small test experiment Mini-Max (T-864) [ES] at Tevatron was specially 
designed for a search of DCC in the forward region. The results of this ex- 
periment j27j are consistent with the generic production mechanism and 
show no evidence of the presence of DCC. Despite this failure to find DCC, 
the Mini-Max experiment was an important benchmark. It was demon- 
strated for the first time that it was possible to work in the very forward 
region with severe background conditions. Much was learned in both detec- 
tor operation and data analysis which should prove useful in future more 
elaborate efforts of this kind. 

Central rapidity region Centauros were searched in the CDF experiment 
at Tevatron with negative result |94j. Another major Tevatron detector DO 
is also suitable for such searches, as the Monte Carlo study shows [94J. 

A serious effort to study possible DCC formation in heavy ion collisions 
was undertaken in the CERN SPS fixed target experiment WA98 |98j. Again 
no DCC signal was found in the central 158 • A GeV Pb+Pb collisions. 

A majority of future heavy ion experiments at RHIC and LHC have plans 
to look for the Centauro phenomena [94] . The kinematic conditions at which 
cosmic ray Centauros are produced will be accessible at RHIC. Therefore 
the corresponding experiments (PHOBOS, STAR, PHENIX and BRAHMS) 
are very interesting in light of Centauro investigation. At LHC the energy 
accessible in Pb+Pb central collisions will be much higher than the expected 
threshold energy for the Centauro production. Besides, Pb+Pb collisions at 
LHC will produce a very high baryon number density in the forward rapidity 
region. To study the novel phenomena expected in such high baryochemical 
potential environment, the CASTOR detector, as the part of the ALICE 
experiment, was designed [99J. Its main goal is the Centauro and strangelet 
search in the very forward rapidity region in nucleus-nucleus collisions. 

We believe that future photon-photon colliders are also good places to 
look for the DCC production. Some hints were given above that the DCC 
formation conditions might be even more favorable at photon-photon collid- 
ers rather than at proton-proton (or proton-antiproton) colliders. Here we 
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present one more argument of this kind which deals with the very different 
roles played by gluons in mesons and baryons. Mesons can be considered 
as a quark- antiquark pair connected by a gluon string (flux tube). There- 
fore, the gluon field configuration in mesons is, in some sense, topologically 
trivial. In baryons one has a quite different picture jlOOl 11011 1102j . Accord- 
ing to the common wisdom, baryons are three-quark bound states. In high 
energy pp or heavy ion collisions the valence quark distributions in the pro- 
jectiles will be significantly Lorentz-contracted because a typical fraction of 
the proton's momentum carried by a valence quark is ~ 1/3. Therefore, 
one expects that the constituent quarks of the colliding protons will not 
have enough time to interact significantly during the high-energy collisions 
and hence it is difficult to stop them. If the baryon number of the projec- 
tile is associated with their valence quarks, which is the naive expectation, 
then the ready prediction from the above given collision picture will be that 
the baryon number flow should be concentrated at large positive and nega- 
tive rapidities, with a nearly zero net baryon number at central rapidities. 
Surprisingly, this is not the case supported by experiments. On the con- 
trary, experiments suggest that the valence quarks do not carry the proton's 
baryon number and the flow of the baryon number can be separated from 
the flow of the valence quarks jlOlL ll()2j ! But then what is the mysterious 
fourth constituent of the proton which traces its baryon number? In QCD 
the baryon is represented by a gauge invariant, non-local, color singlet oper- 
ator. In fact, the gauge invariance constraint severely restricts the possible 
forms of such composite operator, leaving only one possibility (a, (3, 7 are 
the color indices, the flavor indices are suppressed for simplicity) : 



B = e 



T(x h x)q(x 1 ) T(x 2l x)q(x 2 ) T(x 3 ,x)q(x 3 ) 



Here T( 1 . X J IS the open string operator (the Wilson line), or parallel trans- 
porter of the quark field q(xi) from the point X{ to the point x, where the 
three strings join. This string operator is an analog of the well known 
Aharonov-Bohm phase in QED and is given by the path-ordered exponent 




f( Xi , x) = P exp ( ig I A^dx^ ) , A /l = A a ^ 



Therefore, the gluon strings (flux tubes) inside a baryon have nontrivial 
Y-shaped topology and one finds a novel object there - the string junction. 
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This string junction is just the fourth constituent of the baryon which traces 
its baryon number [101J. The string junction can be more easily stopped in 
the high-energy collisions, because, being formed from the soft gluons, it is 
not Lorentz-contracted and always has enough time to interact. 

Now we have the following picture of the high-energy pp collisions jlOlj : 
the valence quarks are stripped-off and produce jets in the fragmentation 
regions. In some events, one or both of the string junctions are stopped 
in the central rapidity region producing a violent gluon sea containing one 
or two twists. On the contrary, no twists are expected in the gluon sea 
produced by high-energy photon-photon (pp) collisions. We believe that 
the latter situation is more favorable for the Baked Alaska scenario and, 
therefore, for the DCC production through this mechanism. 

As a final remark, let us note that the DCC formation is just one in- 
teresting collective effect expected in high-energy collisions. Other exotic 
phenomena are also worth to be searched. Let us mention a few: the possi- 
ble formation of the pion and eta strings during the chiral phase transition 
[103J, creation of the parity and CP violating metastable vacuum bubbles 
[104J, production of QCD Buckyballs - femtometer scale gluon junction net- 
works (QCD analog of the nanoscale carbonic Fullerenes) [102J. Vacuum 
engineering at photon colliders promises to be an exciting adventure and we 
suspect that one may encounter "totally unexpected" new phenomena dur- 
ing such exploration: "There are more things in Heaven and Earth, Horatio, 
than are dreamt of in your philosophy" [105J. 
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Appendix: the Baked Alaska recipe 

Here we reproduce the recipe from [I106| . 

Ingredients: 

• 3 egg whites 
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• 1/2 cup of sugar 

• 1 cup of really hard, frozen ice cream 

• 1 big, thick, hard cookie 

• Baking sheet 

• Aluminum foil 

• Hand mixer 

• A grown up!! 

Directions: 

1. Have your grown up heat your oven to 500 degrees Fahrenheit. 

2. Cover the baking sheet with aluminum foil. 

3. Put the egg whites into a bowl and use the mixer to beat them for about 
five minutes until they're stiff. 

4. Keep beating the egg whites while adding the sugar a little at a time 
until they're fluffy and shiny. ("The name for this stuff is meringue!") 

5. Put your cookie on the baking sheet. 

6. Put your scoop of ice cream on top of the cookie. Make sure it doesn't 
hang over the edge of the cookie. 

7. Completely cover the cookie and the ice cream with the meringue. 

8. Put it in the oven for three to five minutes until the meringue is a delicate, 
light brown. 

9. Take it out of the oven, put it on a plate, and eat up! 
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